Measurement of the interaction strength in a Bose-Fermi mixture with ^''Rb and 
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A quantum degenerate, dilute gas mixture of bosonic and fermionic atoms was produced using 
*''Rb and '^''K. The onset of de generacy was confirmed by observing the spatial distribution of the 
gases after time-of-flight expansion. Further, the magnitude of the interspecies scattering length 
between the doubly spin polarized states of *^Rb and *''K, |aRbK|, was determined from cross- 
dimensional thermal relaxation. The uncertainty in this collision measurement was greatly reduced 
by taking the ratio of interspecies and intraspecies relaxation rates, yielding |aRbK| = 250 ± 30 ao, 
which is a lower value than what was reported in [M. Modugno et al., Phys. Rev. A 68, 043626 
(2003)]. Using the value for |aij,bK| reported here, current T = theory would predict a threshold 
for mechanical instability that is inconsistent with the experimentally observed onset for sudden 
loss of fermions in [G. Modugno et al, Science 297, 2240 (2002)]. 
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With many recent advances, ultracold Fermi gas ex- 
periments are now beginning to explore the strongly in- 
teracting regime, where pairing can lead to fermionic su- 
perfluidity 1]. A natural extension to such studies is 
to investigate the properties of a mixture of a Fermi 
gas and a Bose-Einstein condensate (BEG). The pres- 
ence of a BEG can vastly modify the static and dynamic 
character of a quantum degenerate Fermi gas. M0lmer 
pointed out that the mixture can phase separate, coexist, 
or collapse depending on the relative sign and magnitude 
of the boson-boson and boson-fermion scattering lengths 
. The spectrum of collective excitations can show rich 
physics as the interaction strength between the gases in- 
creases In addition, an effective fermion-fermion at- 
traction arises from the boson-fermion interaction, which 
can significantly increase for a phase transition to a 
BGS-type superfluid 0. This mechanism for Gooper 
pairing of atoms is intriguing because of the analogy to 
superconductivity. 

To date, quantum degenerate atomic Bose-Fermi mix- 
tures have been realized in the following combinations: 
^Li (fermion) with ''Li (boson) ^Li with ^^Na (boson) 
0, and ^'^K (fermion) with «^Rb (boson) 0. With the 
last mixture, pioneering experiments by the LENS group 
revealed that the scattering length between ®^Rb and "^^K 
is large and negative [3,11,^7 10]- At high densities, this 
strong attractive interaction can exceed the Pauli pres- 
sure and cause the whole system to become mechanically 
unstable and collapse The LENS group 

has reported the observation of such a collapse; they see 
a sudden loss of "^''K atoms as the number of atoms in 
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FIG. 1: Schematic of the experimental apparatus. The glass 
cell on the left ("MOT cell") is filled with a background va- 
por of ^'^Rb and ^"K released from alkali-metal dispensers 
attached to the left end of the cell. The cell on the right 
("science cell") is under ultrahigh vacuum, and small coils 
surrounding the cell produce a magnetic trap with strong con- 
finement for evaporative cooling. (The figure does not show 
all the coils.) Atoms are transferred from the MOT cell to 
the science cell using a moving magnetic trap produced by an 
anti-Helmholz coil pair on a translation stage. 



the ^^■^Rb BEG exceeds a threshold value 0. 

Here we report the production of a quantum degener- 
ate Bose-Fermi mixture with *^Rb and ^''K atoms and a 
determination of the magnitude of the cross-species scat- 
tering length. Our experimental apparatus for producing 
the mixture, which requires only a single magneto-optical 
trap (MOT), is outlined in Fig. ^ In brief, our vacuum 
chamber consists of two glass cells — one contains a two- 
species MOT while the other cell (science cell) is used 
for evaporative cooling of the mixture. The lifetime for 
^^Rb atoms in the magnetic trap, which is limited by col- 
lisions with the background gas, is ~ 4 s in the MOT cell 
and more than 100 s in the science cell. The atoms are 
transferred between the cells using a moving magnetic 
trap. 

The experimental procedure for creating the mixture 
is as follows. First, 2 x 10^ ^^Rb atoms and 1 x 10^ "^"K 
atoms are collected in a two-species MOT from the back- 
ground vapor p^ . The atoms are then transferred into 
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a quadrupole magnetic trap, whose magnetic field gradi- 
ent along the symmetry axis in the horizontal direction 
is 99 G/cm. During the last 50 ms of the MOT, the ^'^Rb 
cloud is compressed by reducing the intensity of the hy- 
perfine repumping light. Simultaneously the '^''K cloud is 
Doppler-cooled by reducing the detuning of the trapping 
light to about one natural linewidth (~ 6 MHz). In addi- 
tion, the atoms are optically pumped to their doubly spin 
polarized states, which are |F, Wi?) = |2,2) for ®^Rb and 
|9/2, 9/2} for '**'K. Once in the quadrupole magnetic trap, 
the atoms are transferred to the science cell by physically 
moving the trapping coils, which are mounted on a mo- 
torized translation stage [l^- The observed loss in the 
phase-space density of the ^''Rb cloud during the 81 cm 
(6 sec) travel is only a factor of two. In the science cell, 
the atoms are transferred to an loffe-Pritchard-type mag- 
netic trap for the final cooling. Since there is no MOT 
in the science cell, the cell size can be small (our cell is 
1 cm^ in cross-section) and we can place coils very close 
to the atoms to achieve tight confinement with relatively 
low electric current 27 A). The radial (axial) trapping 
frequency for ^^Rb is Vr,Rh ~ 160 Hz (1^2, Rb = 25 Hz), and 
the trapping frequencies for ''"K are \J rTiRb / tok = 1A7 
times higher, where mRb ("i/f) is the mass of ^^Rb (^°K) 
atoms. 

Simultaneous quantum degeneracy is achieved by radio 
frequency (rf)-induced evaporative cooling of *^Rb and 
sympathetic cooling of *°K. This scheme circumvents the 
problem of vanishing elastic cross section for spin polar- 
ized fermions at ultracold temperatures . The typical 
initial number and temperature of the *^Rb atoms in the 
loffe-Prichard trap before evaporation are A^Rb = 3 x 10* 
and Tpib — 400 /xK; with these initial conditions, we pro- 
duce a nearly pure BEG of 2 x 10''' ^''Rb atoms after 40 
seconds of evaporation. The initial conditions for in 
the loffe-Prichard trap are not known due to the difficulty 
of measuring a small optical density cloud with absorp- 
tion imaging. For sympathetic cooling from 10 fiK to 
0.4 /xK, a range where the cloud can be imaged and 
there is not yet a *^Rb BEG, we observe less than 10% 
number loss with typically 1.4 x 10^ ""^K atoms. 

Quantum degeneracy for the *^Rb gas is evident in the 
formation of a Bose-Einstein condensate, which is charac- 
terized by the emergence of a bimodal distribution in the 
absorption images taken after time-of-fiight expansion. 
The fraction of atoms in the BEG is determined by fit- 
ting the distribution with the sum of the Thomas- Fermi 
profile (inverted parabola) for the condensate and the 
thermal distribution for an ideal Bose gas (Fig. El^A)). 
The temperature is determined from the fit to the non- 
condensed component. For comparison to our data, a 
theory curve that includes both the effects from finite 
number and the positive chemical potential due to boson- 
boson interactions in the condensate is also shown [l7l| . 
The theory and the experiment are in good agreement. 
These data were taken with ^ 1 x 10** ''"K atoms present 
at the end of evaporation. In comparing data taken with 
and without '"^K atoms present, no effect of on the 
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FIG. 2; Quantum degeneracy observed in time-of- flight im- 
ages. (A) Condensate fraction A^bec /iVtotai obtained from 20 
ms time-of-fiight images of *^Rb. The dashed line shows the 
condensate fraction for an ideal Bose gas in the thermody- 
namic limit. The solid line includes the effects from the finite 
number of atoms and the positive chemical potential due to 
boson-boson interactions in the condensate flTII . The phase 
transition occurs at 440 nK. (B) Fugacity of the fermionic 
cloud obtained by fitting the density distribution of ''"K atoms 
after 5 ms time-of-fiight. The fugacity z is equal to the peak 
phase space density of the gas in the classical limit. The dif- 
ferent symbols denote data taken with a *^Rb BEC present 
(filled circles), with a *^Rb thermal cloud present (filled tri- 
angles), and without *^Rb atoms present (open circles). The 
data are systematically shifted from the theory curve expected 
for an ideal Fermi gas (line). A typical value of Tf at T/Tf 
= 0.5 is 220 ± 50 nK. 



* Rb condensate fraction is observed. This is consistent 
with the low "^^K density(^ 3 x 10^^ cm~'^), which results 
in a negligible contribution to the *^Rb mean field. 

The effect of quantum degeneracy on fermions is more 
subtle. The density distribution in time-of-flight becomes 
more "flat-topped" than that for a classical gas due to the 
Pauli pressure. The fugacity z — exp^fix /ksT) (where 
fiK is the chemical potential of ""^K), the cloud size in 
each direction, and the peak density of the gas are used 
as independent parameters for fitting the images. In Fig. 
|2lB) the measured fugacity is plotted as a function of 
T/Tp, where T is determined from the cloud size, and 
Tp is the Fermi temperature. From the fugacity mea- 



3 



surement and assuming an ideal Fermi gas, we find that 
the lowest T/Tp achieved in this experiment was 0.20. 
For comparison to the data, the solid line shows the pre- 
diction for an ideal Fermi gas in a 3D harmonic trap 
p,8]. Most of the data lie to the right of the curve. This 
discrepancy could be explained if our ^''K number cali- 
bration is low by about a factor of two so that T/Tp is 
systematically overestimated by 20%. Note that the fu- 
gacity, which is a measure of the "flatness" , is free from 
calibration errors. The temperature measurement has 
been confirmed by fitting to the thermal wing of ^^Rb 
images taken during the same experimental run. Finally, 
we note that the measured fugacity shows no systematic 
difference between ^"K images taken with and without a 
^^Rb BEC present in the trap. With our highest number 
of 1.8 X 10^ s^Rb atoms in a BEC and 8 x 10"* ''"K atoms 
in the Fermi gas we observe only a slow number loss 1 
sec timescale) and no sudden collapse. 

The strength of inter-species interactions is critical in 
determining the static and dynamic properties of the 
mixture. We have determined the magnitude of the inter- 
species scattering length |aRbK| through measurements 
of cross-dimensional thermal relaxation In general, 
this method has limited accuracy because of systematic 
uncertainties in the atom number, which is typically de- 
termined from absorption imaging. Experimental imper- 
fections such as frequency jitter, imperfect polarization 
of the probe beam, or fluctuation and tilt of the magnetic 
field during imaging tend to reduce probe beam absorp- 
tion. This leads to an underestimate of the number of 
atoms and an overestimate of the cross section. Here 
we avoid this systematic error by taking the ratio of two 
cross section measurements that are both proportional 
to the number of ^^Rb atoms: (1) thermal relaxation of 
'"'K atoms in the presence of *^Rb and (2) relaxation 
of a single-species ^^Rb cloud. We further use the fact 
that the Rb-Rb cross section is already known with high 
accuracy 0. 

The experimental procedure is as follows. First, we 
prepare 1 x 10^ '^''K atoms and 4 x 10^ ^''Rb atoms at 
temperatures T ~ 1/iK in a magnetic trap whose radial 
trapping frequency t'r.Rb is 90 Hz. Both gases are away 
from quantum degeneracy and elastic collisions are in the 
s-wave limit. Multiple 1 ms rf sweeps can be applied to 
reduce the number of *^Rb atoms by as much as a factor 
of four, for the purpose of varying the density. The bias 
magnetic field is then ramped from 4.3 G to 1.4 G to 
increase the radial confinement of the magnetic trap to 
its regular strength {vr^Bh — 160 Hz); this adds energy 
to the radial direction. The duration of the ramp is kept 
long (> 10 ms) compared to the radial trapping period to 
avoid excitation of collective modes, but short compared 
to the relaxation time in the gas. This is always possible 
unless the gases enter into the hydrodynamic regime. 

Just after the ramp, the energy of the gas in the ra- 
dial direction is higher than in the axial direction. The 
cross-dimensional rethermalization time constant is ex- 
tracted by holding the mixture for a variable time in the 
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FIG. 3: Typical relaxation curve for the energy anisotropy. 
Shown is the the square of the ''"K cloud aspect ratio after 5 
ms time-of-flight; this is proportional to the ratio of the ener- 
gies in each direction. The ^"K cloud rethermalizes through 
elastic collisions with *^Rb atoms. The line is a fit to the ra- 
tio of exponentially decaying curves. The error bar represents 
only the statistical uncertainty. 

final trap and measuring the aspect ratio after time-of- 
fiight. Shown in Fig. 13 is a typical relaxation curve for 
the energy anisotropy. Rather than the radial and ax- 
ial sizes, the aspect ratio of the cloud, (Jz/<Jri is used for 
data analysis since it is less sensitive to shot-to-shot vari- 
ations in the absolute temperature of the cloud before the 
compression |2ll |. 

The data analysis is simplified by two facts: (i) due to 
the Fermi statistics, the spin-polarized atoms collide 
only with ^^Rb atoms, and (ii) in our experiment there 
is no net flow of energy between the ^''K and ®^Rb gases. 
Therefore, the thermal relaxation rate of the '"'K cloud, 
l/jK) depends only on the collision rate with ®^Rb atoms 

= ^ ("Rb) O-RbK VYihK . (1) 

Tk P 

Here, (nRb) is the average density of *^Rb atoms in the 
trap, (jRbK = 47rap{^j^ is the elastic cross section for Rb- 
K collisions, and URbK = ySfcijTVTr/i is the thermally 
averaged relative speed between ^^Rb and atoms, 
where ^ = TORb™K/('7iRb + Jtik) is the reduced mass. 
The coefficient of proportionality, [3, can be regarded as 
the number of collisions per ''''K atom needed for thermal 
relaxation. The relaxation rate for a single-species *^Rb 
cloud 1/TB.h is proportional to the Rb-Rb collision rate: 

= — ("Rb) cTRbRb vm^ , (2) 

TRb a 

where fiRbRb — Stto^j^ is the elastic cross section for Rb- 
Rb collisions, VRb = ^i^^JhsT /irmYih is the thermally av- 
eraged relative speed between ®^Rb atoms, and a is the 
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FIG. 4: Measurement of elastic collision cross-section through 
energy anisotropy relaxation. Relaxation time constants for 
a *'^K gas in the presence of ^'^Rb (open circles) and ^'^Rb gas 
alone (filled circles) are shown as a function of the product 
of average *^Rb density in the trap and the square root of 
the equilibrium temperature. Lines are linear fits through 
the origin. The ratio of the slopes is used to extract the 
magnitude of interspecies scattering length |aRbK|, as shown 
in Eq. Q. 

coefficient of proportionality for single-species relaxation. 
Thus, the ratio of the two time constants is free from sys- 
tematic uncertainties in the ^^Rb number p^ . 



gives rise to the difference between a and f3, as a lighter 
particle can redistribute more of its total kinetic energy 
in a single collision. 

The measured relaxation rates for varying density of 
*^Rb atoms are plotted in Fig. 21 The error bars only 
account for statistical errors. Systematic uncertainties in 
the ^^Rb density were eliminated from the determination 
of IflRbKl by taking the ratio of the two slopes and using 
Eq. The known value ORb = 98.98 ± 0.04 ao for ^'^Rb 
in the |2,2) state |23|, obtained by combining multiple 
high-precision measurements including photoassociation 
and Feshbach spectroscopy, is used for calibration [25l | 
and we obtain loRbKl = 250 ± 30ao. 

This measured value of loRbxl is in mild disagreement 
with the value of -410+^° oq from the LENS group ||], 
which was determined from a coUisional measurement 
that had the usual large systematic number uncertainty. 
More serious disagreement is seen when comparing our 
result with the value of —395 ± 15 ao determined from a 
comparison of theory predictions and experimental ob- 
servation of collapse phenomena ^Td\. These results sug- 
gest a need for further investigation of the behavior of 
the mixture close to the collapse. For example the zero 
temperature phase diagram in ref. |l2l | which assumed 
ORbK = — 260ao predicts that a condensate with 2 x 10^ 
atoms remains mechanically stable with numbers of '**'K 
atoms as high as 1 x 10^. The new value for |aRbKl that 
we report here is also likely to have an irnpact on the 
predicted positions of Feshbach resonances |2a, |23| . 



TRb _ 1 a / mRb / flRbK \ 

TK " 2 /3 V 2/^ ^ aRb ; ■ ^ ' 

The constant a {(3) has been determined to be 2.67 ± 
0.13 (2.1 ± 0.2) from Monte Carlo simulations using our 
experimental parameters, including the initial energy 
anisotropy of both species and the mass difference be- 
tween "^^K and ^^Rb |2^. It is this mass difference that 
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